arXiv:1506.06063vl [hep-ph] 19Jun2015 


June 22, 

[JMPA-D-14-00127'arXiv 


2015 0:19 WSPC/INSTRUCTION 


FILE 


International Journal of Modern Physics A 
© World Scientific Publishing Company 


NEUTRINO MIXING WITH NON-ZERO 0 13 AND CP VIOLATION 
IN THE 3-3-1 MODEL BASED ON S 4 FLAVOR SYMMETRY 

VO VAN VIEN 

Department of Physics, Tay Nguyen University, 567 Le Duan, Buon Ma Thuot, Vietnam 

wvienkl 6@gmail. com 

HOANG NGOC LONG 

Institute of Physics, VAST, 10 Dao Tan, Ba Dinh, Hanoi, Vietnam 
hnlong@iop. vast, ac.vn 

DINH PHAN KHOI 

Department of Physics, Vinh University, 182 Le Duan, Vinh City, Nghe An, Vietnam 

khoidp@vinhuni. edu. vn 

Received Day Month Year 
Revised Day Month Year 

The 3-3-1 model proposed in 2011 based on discrete symmetry S 4 responsible for the 
neutrino and quark masses is updated, in which the non-zero $13 is focused. Neutrino 
masses and mixings are consistent with the most recent data on neutrino oscillations 
without perturbation. The new feature is adding a new SU(S)l anti-sextet lying in dou¬ 
blet under S 4 which can result the non-zero $13 without perturbation, and consequently, 
the number of Higgs multiplets required is less than those of other models based on 
non-Abelian discrete symmetries and the 3-3-1 models. The exact tribimaximal form 
obtained with the breaking S 4 —»• Z 3 in charged lepton sector and S 4 —>• JC in neutrino 
sector. If both breakings 54 —» JC and JC —» Z 2 are taken place in neutrino sector, the 
realistic neutrino spectrum is obtained without perturbation. The upper bound on neu¬ 
trino mass and the effective mass governing neutrinoless double beta decay at the tree 
level are presented. The model predicts the Dirac CP violation phase S = 292.45° in the 
normal spectrum (with O 23 / and 5 = 303.14° in the inverted spectrum. 

Keywords : Neutrino mass and mixing; Non-standard-model neutrinos, right-handed neu¬ 
trinos; Flavor symmetries; Discrete symmetries; Models beyond the standard model. 
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1. Introduction 

Nowadays, particle physicists are attracted by two exciting subjects: Higgs and 
neutrino physics. The neutrino mass and mixing are the first evidence of beyond 
Standard Model physics. Many experiments show that neutrinos have tiny masses 
and their mixing is sill mysterious 1 ^ . The tri-bimaximal form for explaining the 
lepton mixing scheme was first proposed by Harrison-Perkins-Scott (HPS), which 
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apart from the phase redefinitions, is given bj@£l 

/ -4 -4 o 

kWs = 


2 1 

VE VE 
_ i _ i_ 

VE VE 
_ 1 _ 1 _ 

VE VE 


( 1 ) 


can be considered as a good approximation for the recent neutrino experimental 
data. 

The most recent data are a clear sign of rather large value 6 *i 3 7 . The data in 
the Particle Data Group PDG201 Z P imply: 


sin 2 ( 26 * 12 ) = 0.846 ± 0.021, sin 2 (2<9 23 ) = 0.999tg;g?g, 
sin 2 (26»i 3 ) = (9.3 ± 0.8) x 10" 2 , A m\ x = (7.53 ± 0.18) x 10" 5 eV 2 , 

Ato 32 = (2.44 ± 0.06) x 10 _ 3 eV 2 , (Normal hierarchy), ( 2 ) 

sin 2 ( 26 * 12 ) = 0.846 ± 0.021, sin 2 (26» 23 ) = 1.000tg;$$, 

sin 2 (26»i 3 ) = (9.3 ± 0.8) x 10~ 2 , A m\ x = (7.53 ± 0.18) x 10~ 5 eV 2 , 

Ato 32 = (2.52 ±0.07) x 10 _ 3 eV 2 , (Inverted hierarchy). (3) 


These large neutrino mixing angles are completely different from the quark mixing 
ones defined by the Cabibbo- Kobayashi-Maskawa (CKM) matrix^®, and they can¬ 
not be explained by the Standard Model. It is an interesting challenge to formulate 
dynamical principles that can lead to the flavor mixing patterns for quarks and lep¬ 
tons given in a completely natural way as first approximations. This has stimulated 
work on flavor symmetries and non-Abelian discrete symmetries are considered to be 
the most attractive candidate to formulate dynamical principles that can lead to the 
flavor mixing patterns for quarks and lepton. There are many recent models based 
on the non-Abelian discrete symmetries, such as A 4 (RefsQTHMI , A 5 (Refs]29H4TJ) 
, 5*3 (Refs|42H83|) , S 4 lRefsl84UTl2ll. D 4 lRefs fTT3HT24l) , D 5 (RefsQl5l[l26]) , T' 
fRefs ll27fil31l) and so forth. In our previous workl l 3 StkM 2 l , the discrete symmetries 
have been explored to the 3-3-1 models. In Ref. 11331 we have studied the 3-3-1 model 
with neutral fermions based on £4 group, in which most of the Higgs multiplets are 
in triplets under £4 except \ lying in a singlet, and the exact tribimaximal form 3 6 
is obtained, where 6 *i 3 = 0 . 

As we know, the recent considerations have implied d 33 76 (jllASMaHnU , but 
small as given in Eqs. 0 and 0 . This problem has been improved in Ref. 11341 by 
adding a new triplet p put in 1 ' under S 3 and another antisextet s 1 put in 2 under 
f> 3 , in which s' is regarded as a small perturbation, or a new triplet p put in 1 " 
under D 4 regarded as a small perturbation 135 . Therefore the models contain up to 
eight Higgs multiplets, and the scalar potential of the model is quite complicated. 

In this paper, we introduce another SU(3)l antisextet lying in 2 under £4 which 
can result the non-zero $i 3 without perturbation. The rest of this work is organized 
as follows. In Sec. [2] we review some main results from Ref. 11331 Sec. [3] is devoted for 
the neutrino mass and mixing. Sec. |4] presents the remark on the vacuum alignments 
and p parameter. We summarize our results in the Sec. [5] |Appcndix A| is devoted 
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to * 5*4 group with its Clebsch-Gordan coefficients. |Appendix B| presents the lepton 
numbers and lepton parities of model particles. |Appendix C| provides the breakings 
of S 4 group by triplets 3 and 3'. 


2. The model 


The fermions in this model under [SU(3 )l, U(1).y, U(l)£, 5 4 ] symmetries, respec¬ 
tively, transform a^ 33 ! 


IpL = ^1,2,3-L = 



l 1R ~ [ 1 ,- 1 , 1 , 1 ] 


Q3L = 



[3,-1/3,2/3,3], 

[ 1 ,- 1 , 1 , 2 ], 

[3,1/3,-1/3,1], 


*1,2,3 R/ 

Ir = h,3 R 


(4) 


Ql = Ql, 2 L = — 


Ur = Wl,2,3fl 



[3*, 0,1/3,2], 


Ar = di,2,3fi ~ [1, —1/3,0,3], 


U R ~ [1,2/3, -1,1], D R = D li2R ~ [1,-1/3,1,2], 

where the numbered subscripts on field indicate respective families and define com¬ 
ponents of their S 4 multiplet representation. Note that the 2 for quarks meets the 
requirement of anomaly cancelation where the last two left-quark families are in 3* 
while the first one as well as the leptons are in 3 under SU(3)l ■ All the C charges 
of the model multiplets are listed in the square brackets. 

To generate masses for the charged leptons, we have introduced two SU(3)l 
scalar triplets (j) and </' lying in 3 and 3' under S 4 , respectively, with the VEVs 
(</>} = (v,v,v) and {</'} = (i/,i/,i/) written as those of S 4 components 1 -^, i.e, S 4 is 
broken into Z 3 that consists of the elements^ {1, T, T 2 } . From the invariant Yukawa 
interactions for the charged leptons, we obtain m e = \f3h\V, = v / 3 (/i 2 V — 

/ 13 ?/), m r = V3(h 2 v + J 13 V'), and the left and right-handed charged leptons mixing 
matrices are given 133 


U L = 



U r = 1 . 


(5) 


In similarity to the charged lepton sector, to generate the quark masses, we 
have additionally introduced three scalar Higgs triplets x, 77 , 7/ lying in X, 3 and 


a With the VEV alignment: = ((f) 2 ) = (4> 3 ) 7 ^ 0, ^4 group is broken into S 3 which consisting of 

the elements {1, T, T 2 , TSTS 2 , STS 2 , S 2 TS}; with the VEV alignment: (</)[ ) = ((j)' 2 ) = (0 ^) / 0, 
54 is broken into Z 3 that consists of the elements {1 ,T, T 2 } as presented in | Appendix C| 
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3' under S 4 , respectively. Quark masses can be derived from the invariant Yukawa 
interactions for quarks with supposing that the VEVs of 77 , rf and \ are it, v! and v x , 
respectively, where u = (77 J), u 1 = ( 1 / 1 ), v x = (xg) and the other VEVs ( 77 °), ( 1 / 3 ), 
and (xi) vanish due to the lepton parity conservation. The exotic quarks get masses 
mjj = f3 v x an d mfli ,2 = f v x ■ The masses of ordinary up-quarks and down-quarks 
are: 

m u = —V3(h u v + h' u v'), m c = —V3{h u v — h' u v'), nit = V3h^u, 

md = V3(h d u + h' d u'), m s = V3(h d u — h' d u'), mb = V3h^v. ( 6 ) 

The unitary matrices, which couple the left-handed up- and down-quarks to those 
in the mass bases, are Uf = 1 and = 1, respectively. Therefore we get the 
quark mixing matrix Uqkm = = 1 . For a detailed study on charged lepton 

and quark mass, the reader is referred to Ref. 11331 In this work, we add a new 
SU(3)l anti-sextet lying in 2 under S 4 responsible for the non- zero 6 13 without 
perturbation which is different from those in Refs. ll33til35l . The vacuum alignments 
and the gauge boson masses and mixings are similar to those in Refs. [13511411 so 
we will not discuss it further in this work. 


3. Neutrino mass and mixing 


In this type of the models, the neutrino masses arise from the couplings of ’tpiipL to 
scalars, where /© transforms as 3* © 6 under SU(3 )l and 1 © 2 © 3 © 3' under 
S 4 . For the known scalar triplets (<^>, (j)', x, ih v')i the available interactions are only 
and , but explicitly suppressed because of the £-symmetry. We 

will therefore propose new SU(3 )l antisextets, lying in either 1, 2, 3, or 3' under 
S 4 , which interact with iplipL to produce masses for the neutrinos. In Ref ] 1331 we 
have introduced two SU(3)l antisextets a,s transform as follows 


a = 


s = 


n 


n \ 

^11 

CT 12 

a 13 

_+ 


_+ 

a 12 

°22 

°23 

n 


n 

a 13 

a 23 

a 33 J 

nO 


„0 \ 

S 11 

S 12 

s 13 \ 


cjH — I - 

„+ 1 

S 12 

S 22 

S 23 

„0 


„0 / 

s 13 

S 23 

s 33 / 


[6*, 2/3,-4/3,1], 


[6*, 2/3,-4/3,3], 


( 7 ) 


with the VEV of s is set as ((si), 0,0) under S 4 , where 

( A s 0 v s 

0 0 0 

v« 0 A s 


( 8 ) 


and the VEV of a is 


/ 0 v„ \ 

(a) = 0 0 0 

\v a 0 A CT / 


( 9 ) 
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With these SU (3 )l anti-sextets, the exact tribimaximal form was obtained, in which 
013 = (P 33 . However, the recent experimental data have implied $13 7 ^ 0 as given 
in Eqs. © and ©■ So that we need to modify the neutrino mass matrix to fit the 
recent data. 

Notice that the VEV alignment as in ©, £4 is broken into a group 
which is isomorphic to Klein four groupPS that consists of the elements /C = 
{1, £ , TSTS 2 , TST}. To obtain a realistic neutrino spectrum, in this work we ad¬ 
ditionally introduce another SU (3 )l anti-sextet (s') which lies in 2 under £4 and re¬ 
sponsible for the breaking K, —> Z 2 . This happens in any case below: (s') = ((s[), 0), 
with 


[KQ<\ 

(si) =000. (10) 

\v' s 0A'J 

The VEV alignment of s' as in m will break /C into Z 2 that consists of the elements 
{1, A 2 } (instead of £4 is broken into another Klein four group^ that consists of the 
elements {1, £ 2 , T£ 2 T 2 , T 2 S 2 T}). 

In calculation, combining both cases we have the Yukawa interactions responsible 
for neutrino mass: 

- P + 7jy(4>L4>L)3S + \z(^ L ) 2 s' + H.C. 

X — — — 

= ^C01Z,^1L + ^2L^2L + ll>3Lll>3L)cr 

+ f [(^2L^3L + ^3L^2l)S! + (V^IL + 4>i L 1p3L)s2 + (^t L 1p2L + $ 2 L^1 l)s 3 \ 

+ | [W’ilV’U, + ^ 2 ^ 2 Li> 2 L + ufy L 1 p 3 L )s' 2 + (tpiL^lL + + UJ 2 ^l^l) s\] 

+ H.c. (11) 

The mass Lagrangian for the neutrinos is given by 

- £™ ass = lx (A a u c 1L u 1L + v a v c 1L N? R + v a N 1R u 1L + A a N 1R N? R 

+ \rjV 2L v 2L + v cr^2L^2R + v cr^2R l '2L + A a N 2R N 2R 
+ A (J p3 L r'3L + v a v 3L N 3R + v a Af 3R iy 3R + A a N 3R N 3R ) 

+ ^ [K{^2L U 3L + V^l v 2l) + V s (^2L^3R + ^3L^2r) 

+ v s (N 2R v 3L + N 3R v 2L ) + A s (N 2R N^ r + N 2R N^ R )] 

+ 2 +v' s vl l^ir + v' s N 1R v 1L + a ’ s N 1r N^ r ) 

+ w(\' s v 2L V 2 l + v' s v 2L N 2R + v' s N 2R v 2 L + A' s N 2R N 2R ) 

+ u 2 (A ' s v 3 i,V 3 L + v' s v 3L N 3R + v' s N 3R V 3 l + A' S N 3R N 3R )] + H.c. (12) 
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We can rewrite the mass Lagrangian for the neutrinos in the matrix form: 


- C ass = \xlMuXL + H.c., 


XL = 


VL 


N' 


R 


M„ = 


M l Ml 
Md Mr 


(13) 


where v = (^i,^ 2 ,^ 3 ) T and N = (Ni, N 2 , Ns) T . The mass matrices are then ob¬ 
tained by 


' aL,D,R + d,L,D,R 0 0 

Ml,D,R= | o 0,L,D,R + UcLl,D,R i>L,D,R 

0 bL,D,R a-L,D,R + <^ 2 dL,D,R / 


where 


a L = \ a x, a D = v a x, a R = A a x , 

bL = A s y, b D = v s y , b R = A s y, (14) 

d L = X' s z, d D = v' s z, d R = A ' s z. 

The VEVs A^g break the 3-3-1 gauge symmetry down to that of the standard 
model, and provide the masses for the neutral fermions N R and the new gauge 
bosons: the neutral Z' and the charged Y ± and X 0,0 *. The \ a ^ s and v„ tS belong to 
the second stage of the symmetry breaking from the standard model down to the 
SU(3)c < 8 > U(1)q symmetry, and contribute the masses to the neutrinos. Hence, to 
keep a consistency we assume that A CTiS v a<s . . The natural smallness of 

the lepton number violating VEVs X a ^ s and iv jS was explained in Ref il331 . Three 
active-neutrinos therefore gain masses via a combination of type I and type II seesaw 
mechanisms derived from USD as 


where 


with 


(A 0 0 \ 

M e s = Ml — MJjM^Md =| 0 Bi D J , 

V 0 D B 2 J 


A — aL + Al 


(a R + do ) 2 
a R + d R 


B i = 


ai + biuj 2 + cite 
a R + d R ~ a R d R - b 2 R 1 


D = 
B 2 


a 2 - b 2 

a R + d R- a R d R ~ b 2 R 1 
a i + biu + cioj 2 

a R + d R ~ a R d R - b R ’ 


(15) 


(16) 


°i = a2 D a R + “^aoidodR — bL)b R ) + a R (b 2 D — dLd R ) — aL(a 2 R — b 2 R + d 2 R ) + aLa R d R , 
b\ = a 2 D d R + a R (d 2 D — dLd R ), 

Ci = 2 d D (a D a R - b D b R ) + (b 2 D + d 2 D )d R - d L (a R -b R + d R ), (17) 

a 2 = OL,b R — 2 bL>(a R aR + d R di j) + a 2 R bL + b R (b 2 D — bLb R + d^ D ) + j, 
b 2 — —a R bodD + aubndL, — anbodn + a R bLd R . 
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We can diagonalize the mass matrix m as follows: 

U^M eS U v = diag(mi, m 2 , m 3 ), 


where 


mi = i (Bi + H 2 + \/4-D 2 + {B 1 - B 2 ) 2 ) 


m 2 = A, 


(18) 


m 3 = i (Bi + S 2 - + (Si - S 2 ) 2 ) , 

and the corresponding eigenstates put in the lepton mixing matrix: 


where 


U u = 


K = 



Bi - B 2 - sjAD 2 + (b 1 -B 2 ) 2 


2D 


The lepton mixing matrix is defined as 

( - 1-A 

u lep = ulu v = 1 


y/3 


Vk^T I 

JK* + 1 

cj(1 — 


1+A 

v'p+l 

7F+T 



(19) 


( 20 ) 


( 21 ) 


It is easily to check that Ul in © is an unitary matrix. So, if U v in (llT)l) is unitary 
then Ui e p in (1211) is unitary. Here, we will only consider real values for I\ since the 
unitary condition of Ui ep . Furthermore, it is worth noting that in the case of the 
subgroup 1C is unbroken, i.e, without contribution of s' (or X' s = v' s = A' s = 0), the 
lepton mixing matrix m being equal to Uhps as given in ©• 

The value of the Jarlskog invariant Jcp , which gives a convention-independent 
measure of CP violation, is defined from ED as 

Jcp = Im[ U 21 U^ 2 U 32 ] = ( 22 ) 

Until now the values of neutrino masses (or the absolute neutrino masses) as 
well as the mass ordering of neutrinos are unknown. The neutrino mass spectrum 
can be the normal hierarchy (|mi| ~ |m 2 | < |m. 3 |), the inverted hierarchy (|m 3 | < 

|mi| ~ |m 2 |) or nearly degenerate (|mi| ~ |m 2 | ~ |m 3 |). An upper bound on the 
absolute value of neutrino mass was found from the analysis of the cosmological 
data(442l 


rm < 0.6 eV, 

while the upper limit on the sum of neutrino masses given in 143 

3 

mi < 0.23 eV 

2—1 


(23) 


(24) 
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In the case of 3-neutrino mixing, the two possible signs of Am | 3 corresponding to 
two types of neutrino mass spectrum can be provided as follows: 

o Normal hierarchy (NH): \m±\ ~ |m 2 | < 1 777-3 1, A m§ 2 = *t*| — m\ > 0. 
o Inverted hierarchy (IH): |m 3 | < \m\\ ~ |m 2 |, Am§ 2 = ***§ — m| < 0. 

As will be discussed below, the model under consideration can provide both normal 
and inverted mass hierarchy. 


3.1. Normal case > 0) 

In the Normal Hierarchy, combining (E^l) with the data in Ref. El J C p = -0.032, 
we get 


K = -1.41297, 


and the lepton mixing matrices are obtained as 


or 


Ulep — 


/ 0.805 A= 0.138 \ 

-0.402 + 0.119* 0.069 + 0.697* 

v -0.402 - 0.119* Ag 0.069 - 0.697 *) 


xP, 


/ 0.805 
\U lep \ = 0.420 

V 0.420 


0.577 0.138 \ 

0.577 0.700 , 

0.577 0.700/ 


(25) 


(26) 


(27) 


In the standard parametrization, the lepton mixing matrix can be parametrized as 


Upmns — 


( C12C13 S12C13 si3e _!<5 \ 

— S 12 C 23 — Ci2S23Si3e i<5 C 12 C 23 — Sl2S23Sl3e l<5 S 23 C 13 I .V, (28) 

S 12 S 23 — Ci2C23Sl3e l<S —C 12 S 23 — Si2C23S\3e lS C 23 C 13 / 


where V = diag(l, e*“, e i/3 ), and c^- = cos%, Sy = sinOij with (?i 2 , O 23 and ^13 
being the solar, atmospheric and reactor angles, respectively. <5 = [0,27r] is the 
Dirac CP violation phase while a and /3 are two Majorana CP violation phases. 
Using the parametrization in Eq. (l28l) we get 


(29) 


With the help of (P5l) and (E!H) we have sin 5cp = —0.9242, i.e, Scp = —67.55° 
or 5 C p = 292.45°. 

From Eqs. m and (ESI) we get 


Bi = B 2 - 0.705241D. (30) 

In the normal case, i.e, Am § 2 = to§ — m| > 0, taking the central values of neutrino 
mass squared difference from the data in[ 8 ]as shown in EJ : Am|i = 7.53 x 10 -5 eV 2 
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and Am§ 2 = 2.44 x 10 3 eV 2 , with 7711 , 2,3 given in Eq. m , we get a solutiorEI (in 
[eV]) 


B 2 = -0.5a/4 A 2 - 0.0003 - 0.707729D, 

D = 0.471543 (\/H 2 + 2.44 x 10 ” 3 - \] A 2 - 7.53 x 10“ 5 ) . (31) 

With Bi y2 and D in Eqs. (1311) and (l30ll . 7711 , 2,3 depends only on one parameter A , so 
we will consider mi, 2,3 as functions of A. By using the upper bound on the absolute 
value of neutrino mass in (l23l) we can restrict the values of A: | A\ < 0.6 eV. However, 
in this case, A G (0.0087,0.05) eV or A G (—0.05,—0.0087) eV are good regions of 
A that can reach the realistic neutrino mass hierarchy. 

In Fig. [I] we have plotted the absolute value |mi, 2 , 3 | as functions of A with 
A G (0.0087,0.05) eV and A G (—0.05,—0.0087) eV, respectively. This figure shows 
that there exist allowed regions for values of A where either normal or quasi¬ 
degenerate neutrino masses spectrum is achieved. The quasi-degenerate mass hi¬ 
erarchy is obtained if |H| G [0.05 eV, + 00 ). However, |H| must be small enough 
because of the scale of |mi, 2 , 3 |). The normal mass hierarchy will be obtained if A 
takes the values around (0.0087, 0.05) eV or (—0.05, —0.0087) eV. The sum of neu¬ 
trino masses in the normal case = Si=i l TO *l with A G (0.0087, 0.05) eV is 
depicted in Fig. [2] which is consistent with the upper limit given in Eg. (1241) . 


H.WiWCO H>NW( e *0 



Fig. 1. |mi } 2 , 3 | as functions of A in the case of Am^ > 0 with a) A E (0.00867, 0.05) eV and b) 

A e (-0.05, -0.00867) eV. 


b In fact, this system of equations has four solutions, however, these equations differ only by the 
sign of mi ? 2,3 that it is not appear in the neutrino oscillation experiments. So, here we only consider 
in detail the solution in ED- 
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iV) 



Fig. 2. as a function of A with A E (0.00867, 0.05) eV in the case of Am^ > 0. 


From the expressions ([20 |) , ET[) , ([30]) and (f3l]) , it is easily to obtain the effective 
masses governing neutrinoless double beta decajU^M^ , 

m ee=Yl U ei Kl > m P = ( 52 \ U ^ TO i 

i=l \i=l 

which is plotted in Fig. [3] with A € (0.0087,0.05) eV in the case of Ato| 2 > 0. We 
also note that in the normal spectrum, |mi| « |rn 2 | < |to 3 |, so toi given in (fT51) is 
the lightest neutrino mass, which is denoted as m i = m^ ght . 



UK. 


’\ m /3 \’\ rn light 


,JV 


N 


IK 


(er) 


Fig. 3. \m£,\, 

> 0. 



To get explicit values of the model parameters, we set ^4 = 10 2 eV, which is 
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safely small. The other physical neutrino masses are explicitly given as 

|mi| ~ 4.97 x 10~ 3 eV, \m 2 \ = 10" 2 eV, |m 3 | ~ 5.04 x 10“ 2 eV. (33) 


It follows that 


m* | ~ 7.50 x 10“° eV, |m J s v | = 9.87 x 10“° eV 


1-3 


N | 


1 - 3 , 


(34) 

B x = -3.523 x 10“ 2 eV, B 2 = -2.013 x 10“ 2 eV, D = 2.142 x 10“ 2 eV. (35) 

This solution means a normal mass spectrum as mentioned above. Furthermore, by 
assuming thalB 

A s = A' = A ct = 1 eV, v s = v' s = Vo , A’ s = A a = A s , A s = an 2 , 


we obtain a solution 


x ~ (2.0 + 0.2*) x 10~ 3 , y ~ -(6.1 + 0.61*) x 10“ 3 , 
z = -(4.85 + 0.48) x 10“ 3 , a ~ 0.222 + 0.017*. 


(36) 


(37) 


3.2. Inverted case (Am^ < 0^ 

For inverted hierarchy, the data in Ref. [§] implies Jcp = —0.029. Hence, we get 

K = -1.36483, (38) 


and the lepton mixing matrices are obtained as 


or 


Ui 


ep — 


/ 0.807 Aj -0.125 \ 

-0.403 + 0.108* Aj 0.062 + 0.699* 
v -0.403 - 0.108* Aj 0.062 - 0.699 *) 


xP, 


/ 0.807 
\U lep \ = 0.418 

V 0.418 


0.577 0.125 \ 

0.577 0.701 

0.577 0.701 / 


Combining | 5 ]), ( 1551 ) and ( 1551 ) yields sin Jcp = —0.8371, i.e, Scp 
5 CP = 303.14°. 


From Eqs. m and ( 1551 ) we get 


(39) 


(40) 
-56.84° or 


Bx = B 2 - 0.632138 D. (41) 

In the inverted case, Am§ 2 = to 2 —m 2 < 0, taking the central values of neutrino mass 
squared difference from the data in Ref.IS]as shown in ([5]): Am^ = 7.53 x 10“ 5 eV 2 


c The values of the parameters A s , A^, v s , v ' s , A s , A' s , have not been confirmed by exper¬ 
iment, however, their hierarchies were given in Ref. 11341 The parametres in Eqs. (I36l and ED is 
a set of the model parameters that can fit the experimental data on neutrino given in ©• 
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and A m§ 2 = —2.52 x 10~ 3 eV 2 , with 7711 , 2,3 given in Eq. m, we get a solutiorU 
(in [eV]) 

B 2 = —0.5\/4A 2 - 0.0003 - 0.732692 D, 

D = -0.476753 (yA 2 + 2.52 x 10" 3 + \J A 2 - 7.53 x 10“ 5 ) . (42) 

With £>i ,2 and D in Eqs. m and (1421) . toi, 2,3 depends only on one parameter 
A, so we will consider 7771,2,3 as functions of A. In this case, A £ (0.05, 0.1) eV or 
A £ (—0.1,—0.05) eV are good regions of A that can reach the realistic neutrino 
mass hierarchy. 

In Fig. U we have plotted the absolute value |rrti, 2 , 3 | as functions of A with 
A £ (0.05,0.1) eV and A £ (—0.1,—0.05) eV, respectively. We see that there exist 
allowed regions for values of A where either inverted or quasi-degenerate neutrino 
masses spectrum is achieved. The quasi-degenerate mass hierarchy is obtained if 
|A| £ [0.1 eV,+ 00 ). However, |A| must be small enough because of the scale of 
|t 77 i, 2 , 3 1. The inverted mass hierarchy will be obtained if |A| takes the values around 
(0.05, 0.1) eV. The sum of neutrino masses in the normal case = Y^=i l TO *l with 
A £ (0.0087, 0.05) eV is depicted in Fig. [2] which is consistent with the upper limit 
given in Ea. (l24ll . The effective masses governing neutrinoless double beta decay 


| m 1 |>| m 2 |>| m 3|( e O H|>| m 2|’| ffl 3|( e O 



A e (-0.1,-0.05) eV. 


defined in (l32l) is plotted in Fig. [6] with A £ (0.05,0.1) eV in the case of Am§ 2 < 0. 
We also note that in the inverted spectrum, |to 3 | ~ 1 777 . 2 1 — |tt7i|, so m 3 given in 
m is the lightest neutrino mass, which is denoted as m 3 = mf^ t . 

d Similarly to the normal case, there are four solutions in the inverted hierarchy. Here we only 
consider in detail the solution in ll42)l . 


















June 

[JMPA- 


22, 2015 0:19 WSPC/INSTRUCTION FILE 

i-14-00127'arXiv 


XV) 



Fig. 5. as a function of A with A E (0.05, 0.1) eV in the case of < 0. 


m llght - m 'ee , m 'p ^V) 



Fig. 6. nil and as functions of A with A E (0.05, 0.1) eV in the case of Amg 2 < 0. 


To get explicit values of the model parameters, we set A = 5 x 10 2 eV, which 
is safely small. The other physical neutrino masses are explicitly given as 

|toi| ~ 4.924 x 10 _2 eV, \m 2 \ = 5 x 10" 2 eV, \m 3 \ ~ 4.472 x 10" 3 eV. (43) 

It follows that 

\m{ e \ ~ 4.88 x 10 -2 eV, |m£| = 4.91 x 10" 2 eV, (44) 

Si = (1.72 - 0.29*) x 10“ 2 eV, B 2 = (3.204 - 0.156*) x 10” 2 eV, 

D = (2.348 + 0.213*) x 10“ 2 eV. (45) 
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This solution means an inverted mass spectrum. Furthermore, by assuming that |j 
As — A s = Xa = o, Vg - v s — U(j, A s = A s — A ct , 

As = v 2 s , A' = v' s 2 , A,, = -vl, (46) 


we obtain a solution 

x ~ (3.192 - 0.452*) x 10" 2 , y ~ (-2.563 + 0.294i) x 10“ 2 , 
z = -(1.910 + 0.825*) x 10“ 2 , a ~ 0.105 - 0.186*. (47) 


4. Remark on the vacuum alignments and p parameter 


In the model under consideration, to generate masses for all fermions, we need 
eight Higgs scalars 0, 0 ',x,? 7 , ? i',<j,s,s'. It is important to note that \ and s' do 
not break £>4 since they are put in 1 under £>4 while s ', 0 , rj; 4 >', rf can break £>4 into 
its subgroups since they are put in non-trivial representations 2,3,3' of £ 54 . The 
breaking of £54 group depends on the vacuum alignment of the flavons. 

For doublets 2 (s') we have two followings alignments. The first alignment, 0 4 
(s[) 4 (si j) = 0 or 0 / (s' 2 ) 4 (s' 4 ) = 0 or 0 / (s',) 4 ( s ' 2 ) 4 0 then £>4 is broken 
into a group which is isomorphic to Klein four group 7 that consists of the elements 
{ 1 , TS 2 T 2 , £> 2 , T 2 S 2 T}. The second alignment, (s,) = (si )) 4 0 then £>4 is broken 
into D 4 consists of the elements {1, TSTS 2 , TST, S, S 3 , TS 2 T 2 1 S 2 , T 2 S 2 T}. For 
triplets 3 and 3' the breakings of £>4 are given in |Appcndix C[ 

To obtain a realistic neutrino spectrum, in this work, we argue that the breaking 
S4 —1 Z 3 is taken place in charged lepton sector while both breakings £54 — > K, and 
/C —y Z 2 must be taken place in neutrino sector. 

Note that A CT , A s , A CT ' are needed to the same order and not to be so large that 
can naturally be taken at TeV scale as the VEV v x of x- This is because v a ,v s 
and v a > carry lepton number, simultaneously breaking the lepton parity which is 
naturally constrained to be much smaller than the electroweak scale 132 134 150 151 
. This is also behind a theoretical fact that v x , Act, A s , A CT ' are scales for the gauge 
symmetry breaking in the first stage from SU(3 )l (g) U(1 ).y —> SU(2)l (g U(l)y in 
the original form of 3-3-1 models 150 153 . They provide masses for the new gauge 
bosons Z\ X and Y. Also, the exotic quarks gain masses from v x while the neutral 
fermions masses arise from A CT ,A S ,A CT '. The second stage of the gauge symmetry 
breaking from SU(2)i®U(l)y —> U(1)q is achieved by the electroweak scale VEVs 
such as u,v responsible for ordinary quark masses. In combination with those of 
type II seesaw as determined, in this type of the model, the following limit is often 

taken into accounP^HIMltllMlII . 


(eV ) 2 


\2 \2 ,, 2 2 2^2 /2 2 !\ 
Act> A s , Act/ < v a ,v s ,v a , < u ,u ,v ,v 


< v( 


A 2 ~ A 2 ~ A 2 , ~ (TeV) 2 . 


(48) 


e The values of the parameters A s , A', Ao-, v s , v ' s , v G , A s , , A CT have not been confirmed by exper¬ 
iment, however, their hierarchies were given in Ref. 11341 The parametres in Eqs. (QU and G3 is 
a set of the model parameters that can fit the experimenta data on neutrino given in ®. 
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On the other hand, our model can modify the precision electroweak parameter 
such as p parameter at the tree-level. To see this let us approximate the masses of 
W and Z bosons 0: 


M w - 2g 2 (3u 2 - i%), 

My - y ( QA l + 4A "' + 2 A l' + v l) ■ 

The p parameter is defined as 


P = 


Mh 
c 2 w M 2 z 



(49) 

(50) 


(51) 


It is easily to see that the p parameter in (ICT1) is absolutely close to the unity since 
v 2 <C u 2 and this is in agreement with the data in Ref{Sj 

The mixings between the charged gauge bosons W — Y and the neutral ones 
Z' — W4 are in the same order since they are proportional to j 2 -, and in the limit 
v a <C Act these mixing angles tend to zero. In addition, from C51) and (Hill) . (1501) . it 
follows that M^r is much smaller than My. 


5. Conclusions 

In this paper, we have modified the previous 3-3-1 model combined with discrete 
£4 symmetry to adapt the most recent neutrino mixing with non-zero $ 13 . We have 
shown that the realistic neutrino masses and mixings can be obtained if the two di¬ 
rections of the breakings £4 —> K, and K, —> Z 2 simultaneously take place in neutrino 
sector and are equivalent in size, i.e, the contributions due to s, er and s' are compa¬ 
rable. The new feature is adding a new SU(3)l anti-sextet lying in 2 under £>4 which 
can result the non-zero 9 \3 without perturbation, and consequently, the number of 
Higgs multiplets required is less than those of other models based on non-Abelian 
discrete symmetries and the 3-3-1 models. The exact tribimaximal form obtained 
with the breaking £4 —> Z 3 in charged lepton sector while £4 —> K, in neutrino 
sector. If both the breakings £4 —> K, and K, —> Z 2 are taken place in neutrino 
sector, the realistic neutrino spectrum is obtained without perturbation. The upper 
bound on neutrino mass as well as the effective mass governing neutrinoless double 
beta decay at the level are presented. The model predicts the Dirac CP violation 
phase <5 = 292.45° in the normal spectrum (with $23 7 ^ f) and S = 303.14° in the 
inverted spectrum. We have found some regions of model parameters that can fit 
the experimental data in 2014 on neutrino masses and mixing without perturbation. 


f We have used the notation % = sin (i\\- , c\y = cos 0 w , t\y = tan , and the continuation of the 
gauge coupling constant g of the SU(3) c at the spontaneous symmetry breaking poi ni l 1 3s| EI 3 1 1 1 v J 'l 
3y/2s\y 


t = 




was used. 
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Appendix A. Sa group and Clebsch-Gordan coefficients 

£4 is the permutation group of four objects, which is also the symmetry group of a 
cube. It has 24 elements divided into 5 conjugacy classes, with 1, 1', 2, 3, and 3' as its 
5 irreducible representations. Any element of S 4 can be formed by multiplication of 
the generators S and T obeying the relations S A = T 3 = 1, ST 2 S = T. Without loss 
of generality, we could choose S = (1234), T = (123) where the cycle (1234) denotes 
the permutation (1,2,3,4) —> (2,3,4,1), and (123) means (1,2, 3,4) —► (2, 3,1,4). 
The conjugacy classes generated from S and T are 

Ci : 1 

C 2 : (12) (34) = TS 2 T 2 , (13)(24) = S' 2 , (14)(23) = T 2 S 2 T 
C 3 : (123) = T, (132) = T 2 , (124) = T 2 S 2 , (142) = S 2 T, 

(134) = S 2 TS 2 , (143) = STS, (234) = S 2 T 2 , (243) = TS 2 
C A : (1234) = S, (1243) = T 2 ST, (1324) = ST, 

(1342) = TS, (1423) = TST 2 , (1432) = S 3 
C 5 : (12) = STS 2 , (13) = TSTS 2 , (14) = ST 2 , 

(23) = S 2 TS, (24) = TST, (34) = T 2 S 


The character table of S4 is given as follows 


Class 

n 

h 

Xl 

Xi' 

X2 

X3 

X3' 

Ci 

1 

1 

1 

1 

2 

3 

3 

c 2 

3 

2 

1 

1 

2 

-1 

-1 

c 3 

8 

3 

1 

1 

-1 

0 

0 

Ci 

6 

4 

1 

-1 

0 

-1 

1 

C 5 

6 

2 

1 

-1 

0 

1 

-1 


where n is the order of class and h is the order of elements within each class. 
Let us note that ( 7 ^ 2,3 are even permutations, while C' 4.5 are odd permutations. 
The two three-dimensional representations differ only in the signs of their C 4 and 
C 5 matrices. Similarly, the two one-dinrensional representations behave the same. 
We will work in the basis where 3,3' are real representations whereas 2 is com- 




June 

IJMPA- 


22, 2015 0:19 WSPC/INSTRUCTION FILE 

i-14-00127'arXiv 


17 


plex. One possible choice of generators is given as follows 


1:5=1, T = 1 
l':S=-l, T= 1 



/-I0 0 \ /0 01\ 

3:5= 0 0-1, T = I 1 0 0 J 

\ 0 1 0 / \0 1 0 / 


3 ' :S=- 



T = 



where u> = e 27r? / 3 = —1/2 + *\/3/2 is the cube root of unity. Using them we calculate 
the Clebsch-Gordan coefficients for all the tensor products as given below. 

First, let us put 3(1,2,3) which means some 3 multiplets such as x = 
(xi,X 2 ,X 3 ) ~ 3 or y = (yi, j/ 2 , 2 / 3 ) ~ 3 or so on, and similarly for the other repre¬ 
sentations. Moreover, the numbered multiplets such as (..., ij,...) mean (,...) 
where Xi and yj are the multiplet components of different representations x and 
y, respectively. In the following the components of representations in l.h.s will be 
omitted and should be understood, but they always exist in order in the components 
of decompositions in r.h.s: 


101 = 1(11), ror = 1(11), lor=r(ii), 

102 = 2 ( 11 , 12 ), 1 ' 02 = 2 ( 11 ,- 12 ), 

10 3 = 3(11,12,13), 1' O 3 = 3'(11,12,13), 

10 3' = 3'(11,12,13), l' 0 3' = 3(11,12,13), 

2 0 2 = 1(12 + 21 ) © 1/(12 - 21 ) © 2 ( 22 , 11 ), 

2 O 3 = 3 ((1 + 2)1, w(l + w2)2, w 2 (l + w 2 2)3) 

©3' ((1 - 2)1,07(1 - w2)2, w 2 (l - w 2 2)3) 

2 O 3' = 3' ((1 + 2)l,o;(l + w2)2, w 2 (l + w 2 2)3) 

©3((1 - 2)1, 07(1 - w2)2,w 2 (1 - w 2 2)3) , 

3 0 3 = 1(11 + 22 + 33) © 2(11 + w 2 22 + w33,11 + w22 + w 2 33) 

©3 S (23 + 32,31 + 13,12 + 21) © 3' a (23 - 32,31 - 13,12 - 21), 
3'03' = 1(11 + 22 + 33) © 2(11 + w 2 22 + u 33,11 + w22 + w 2 33) 

©3 S (23 + 32,31 + 13,12 + 21) © 3' a (23 - 32,31 - 13,12 - 21), 
3 0 3' = 1/(11 + 22 + 33) © 2(11 + w 2 22 + w33, -11 - w22 - w 2 33) 

©3(,(23 + 32,31 + 13,12 + 21) © 3 a (23 - 32,31 - 13,12 - 21), 


where the subscripts s and a respectively refer to their symmetric and antisymmetric 
product combinations as explicitly pointed out. We also notice that many group 
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multiplication rules above have similar forms as those of S 3 and A 4 groups. 

In the text we usually use the following notations, for example, (xy')^ = \xy '}3 = 
{x 2 V 3 — x 3 y 2 i x 3 y'i ~ x iy 3 , x iy 2 — x 2 y'l) which is the Clebsch-Gordan coefficients of 
3 a in the decomposition of 3 ( 8 > 3', where as mentioned x = (xi,x 2 ,X3) ~ 3 and 
y’ = {y'l,V2,y'3) ~ 3 '. 

The rules to conjugate the representations 1, V, 2, 3, and 3' are given by 

2*(1*,2*) = 2(2*, 1*), 1*(1*) = 1(1*), r*(l*) = l'(l*), 

3*(1*,2*,3*) = 3(1*, 2*, 3*), 3'*(1*, 2*, 3*) = 3'(1*, 2*, 3*), 

where, for example, 2*(1*,2*) denotes some 2* multiplet of the form (a;*,^) ~ 2*. 

Appendix B. The numbers 

In the following we will explicitly point out the lepton number (L) and lepton parity 
(Pi) of the model particles (notice that the family indices are suppressed): 



Particles 

L 

Pi 

Nr, 

U , d, ^,$'4, <f> 2 ,fa 2 , Vlrfl, % Xt a 33 ’ s 33 

0 

1 

VL, 1 , 

U, D*, ^,<^3 + , rii,V3°, x°l*, xt, 0 - 13 . ° 23 > s 13 > ^23 

-1 

-1 


(T° rr + <?++ 

u ll? u 12 > u 22 5 * 11 ? * 12 ’ *22 

-2 

1 


Appendix C. The breakings of S 4 by triplets 3 and 3 7 

For triplets 3 we have the followings alignments: 

(1) The first alignment: (fa) 7 ^ (fa) 7 ^ (fa) then S 4 is broken into {1} = {identity}, 
i.e. S 4 is completely broken. 

(2) The second alignment: 0 ^ (i)>i) / ( 1 ) 2 ) = (<^ 3 ) ^ 0 or 0 7 ^ (fa) = (fa) 7 ^ (fa) 7 ^ 
0 or 0 ^ (<f> 1 ) = (fa) 7 ^ (<^ 3 } 7 ^ 0 then S 4 is broken into Z 2 which consisting of 
the elements {1,TSTS 2 } or {1,TS'>S' 2 } or {1,S 2 TS}, respectively. 

(3) The third alignment: (fa) = (</> 2 ) = (fa) 7 ^ 0 then S 4 is broken into S 3 which 
consisting of the elements {1, T, T 2 ,TSTS 2 , STS 2 , S 2 TS}. 

(4) The fourth alignment: 0 = (fa) 7 ^ (fa) = (fa) 7 ^ 0 or 0 = (fa) 7 ^ (fa) = (fa) ^ 
0 or 0 = (cj> 3 ) 7 ^ (fa) = (<j>2) 7 ^ 0 then S4 is broken into Z 2 which consisting of 
the elements { 1 ,TSTS 2 } or (1, X 1 5'5' 2 } or {1,S' 2 TS’}, respectively. 

(5) The fifth alignment: 0 = ((f> 2 ) ± (fa) 7 ^ (fa) 7 ^ 0 or 0 = (fa) 7 ^ (fa) 7 ^ (fa) ± 0 
or 0 7 ^ (fa) 7 ^ (fa) 7 ^ (fa) = 0 then £4 is completely broken. 

( 6 ) The sixth alignment: 0 7 ^ (fa) ± (fa) = (fa) = 0 or 0 7 ^ (fa) 7 ^ (fa) = (fa) = 0 
or 0 7 ^ (fa) 7 ^ (fa) = (fa) = 0 then S 4 is broken into Klein four group K, which 
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consisting of the elements {1, S 2 , TSTS 2 , TST} or {1, TS 2 T 2 , STS 2 , T 2 S} or 
{1, T 2 S 2 T, ST 2 , S 2 TS}, respectively. 

For triplets 3' we have the followings alignments: 

(1) The first alignment: (cp}) ^ (cp 2 ) 7 ^ (cp 3 ) then S 4 is broken into {1} = {identity}, 
i.e. 54 is completely broken. 

(2) The second alignment: 0 7 ^ (cp}) / (c p 2 ) = (<t>' 3 ) 7 ^ 0 or 0 7 ^ (cp}) = {<p’ 3 ) 7 ^ (<p 2 ) ± 
0 or 0 / (fp'i) = {<p 2 ) / {(p 3 ) 7 ^ 0 then S 4 is broken into { 1 } = {identity}, i.e. S 4 
is completely broken. 

(3) The third alignment: (cp}) = (cp' 2 ) = (cp 3 } 7 ^ 0 then S 4 is broken into Z 3 that 
consists of the elements {1 ,T,T 2 }. 

(4) The fourth alignment: 0 = (cp' 2 ) 7 ^ (cp}) = (cp' 3 ) 7 ^ 0 or 0 = (cp}) 7 ^ (cp’ 2 ) = (cp' 3 ) ^ 
0 or 0 = (cp 3 } 7 ^ {cp}) = (cp 2 } 7 ^ 0 then S 4 is broken into Z 2 which consisting of 
the elements {I,! 12 !?} or {l,TST} or {1,ST 2 }, respectively. 

(5) The fifth alignment: 0 = {<j)' 2 ) 7 ^ {cp}) ^ {<p' 3 ) ± 0 or 0 = (<p{) 7 ^ (<p' 2 ) 7 ^ {(p' 3 ) ± 0 
or 0 7 ^ (t p}) 7 ^ (< p 2 ) 7 ^ (< f> 3 ) = 0 then S 4 is completely broken. 

( 6 ) The sixth alignment: 0 7 ^ (cp}) ± (cp' 2 ) = (cp' 3 ) = 0 or 0 ^ (cp' 2 ) 7 ^ {(p 3 ) = {cp}) = 0 
or 0 7^ (cp 3 ) 7^ (< p }) = {< p }) = 0 then £>4 is broken into a four-element subgroup 
generated by a four-cycle, which consisting of the elements {1, S', S' 2 , S' 3 } or 
{1,TST 2 ,ST,TS 2 T 2 } or {1, TS, T 2 ST, T 2 S 2 T}, respectively. 
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